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Summary 



O; 

, For normal canonical models, and more generally a vast array of general spherically symmetric location- 

h ■ scale models with a residual vector, we consider estimating the (univariate) location parameter when it is 

I—} ■ lower bounded. We provide conditions for estimators to dominate the benchmark minimax MRE estimator, 

\ and thus be minimax under scale invariant loss. These minimax estimators include the generalized Bayes 

CN 1 estimator with respect to the truncation of the common non-informative prior onto the restricted parameter 

■ space for normal models under general convex symmetric loss, as well as non-normal models under scale 
f-H . invariant LP loss with p > 0. We cover many other situations when the loss is asymmetric, and where other 

\ generalized Bayes estimators, obtained with different powers of the scale parameter in the prior measure, 

" are proven to be minimax. We rely on various novel representations, sharp sign change analyses, as well 

C$ . as capitalize on Kubokawa's integral expression for risk difference technique. Several other analytical 

£h I properties are obtained, including a robustness property of the generalized Bayes estimators above when 

l— ~~ '■ the loss is either scale invariant L p or asymmetrized versions. Applications include inference in two-sample 

. normal model with order constraints on the means. 
>" 
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2i 1. Introduction 



Preamble 



We begin with the normal model in canonical form 

X ~ N(fi,a 2 ),S 2 ~ a 2 xl, independent (n > 1), (1) 

which plays a central role in both statistical theory and practice. Consider situations where addi- 
tional information on (/i, a) is available in terms of parametric restrictions. Bayesian inference in 
such restricted parameter space problems does not, conceptually, present any difficulties as both 
the prior and the resulting posterior will be adapted and will adapt to the constraints. Assess- 
ing the frequentist performance of Bayesian estimators in such situations is, however, considerably 
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more challenging. Such assessments may include, for instance, testing for minimaxity, an evalua- 
tion in comparison to a benchmark procedure such as minimum risk equivariant (MRE) estimator 
or a maximum likelihood estimator (mle), or a study of the frequentist performance of associated 
Bayesian confidence intervals. 

As an illustration, consider model ([1]) with known a and the nonnegative mean restriction /x > 0. 
Despite early discoveries by Katz (1961) and Sacks (1963) that the generalized Bayes estimator 
with respect to the flat prior on [0, oc) is minimax and dominates the MRE estimator Sq(X, S) = X 
under squared error loss, despite various generalizations to other models and location invariant 
losses (Farrell, 1964; Kubokawa, 2004; Marchand and Strawderman, 2005), no other Bayes minimax 
estimators were known until the Maruyama and Iwasaki (2005) findings which provide other Bayes 
minimax estimators under squared error loss. Even then, little has been obtained for estimating 
H in ([1]) for /i > and unknown o. In this case, Kubokawa (2004) obtained, for scale invariant 
squared error loss, a class of minimax improvements on <5 , which includes the generalized Bayes 
estimator S no (X, S) with respect the truncation of the usual non-informative prior onto the restricted 
parameter space (see expression EJ) . 

Our main motivation for his work has been to generalize and better understand Kubokawa's findings. 
The paper consists of various extensions with respect to the loss, the model, and the prior; which 
bypass in a unified way the specific normal case-squared error loss calculations by Kubokawa. 
Several new technical aspects have been developed to meet such challenges. 



1.2. The problem 

As an extension of model (jTJ, we consider spherically symmetric models for an observable (X, U) = 
(X, Ui, . . . , U n ) with density proportional to 

1 fr {x- fi) 2 + \\u\\\ 

and with n > 1, \i > 0, a > 0. The function / : M + — > M + is known, and it assumed throughout 
that: 

tf'Ct) 

f < 0, and —, . decreases in t for t > . (3) 
f(t) 1 ; 

Hereafter, for conciseness, reference to model (j2J) shall be understood to encompass these assump- 
tions on /. Multivariate (for X) versions of (j2J) have been previously considered, namely in recent 
work where robust minimax generalized Bayes estimators of \i without constraints are provided 
(see Fourdrinier and Strawderman, 2010). Various other features of such models are described in 
Section 2.1. 

We consider estimating \i where it is assumed that (//, a) G 6 = {(/x, a) : /i > 0,<r > 0} under 
location and scale invariant loss 

P(^), (4) 
o 

with (i) p absolutely continuous a.e., (ii) p strictly bowled shaped with pit) > p(0) = for all 
t G K, p' < on (—oo,0) and p' > on (0, oo). We also assume that the pair (/, p) leads 
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to risk finiteness, namely that there exists a unique minimum risk equivariant estimator for the 
unconstrained problem. In such cases, it is given by 5o(X,S) = X + cqS with constant risk 
R((p, a), 5q) = E 0) i[p(X + CqS)], and with (also see Remark [3]) 

c = argmin c {£ 0i i[p(X + cS)]} , (5) 

which is uniquely determined by Eq ) i[S(/(X + coS)} = . It is also worth pointing out that Co = 
for symmetric losses p, and consequently that the MRE estimator coincides with the unbiased 
estimator X, and is robust with respect to the choice of the underlying model density /. It follows 
from Kiefer (1957) that Sq is minimax for the unconstrained problem. With the constraint on 
p, 5q(X,S) produces indeed implausible estimates, but it remains minimax (see Marchand and 
Strawderman, 2011, and references therein) for general p, and its constant risk thus matches the 
minimax risk. The challenge here is to search for good improvements on 5q(X,S) that capitalize 
on the parametric information, and we focus on potential Bayesian improvements such as the 
generalized Bayes estimators 5 n with respect to the prior measures 

^O"' a ) = 1 [o,oo)(/^)l(o,oo)(o-) ; I > -(n - 1) ; (6) 

the lower bound on / required for the posterior density to be well defined. The class includes the 
choice 7r which is of intrinsic interest as it represents a plausible adaptation, or truncation onto 
of the right Haar invariant measure 7r r /j with the MRE estimator (also) being the generalized Bayes 
estimator 5 7Trh with respect to ir rh . Moreover, the study of frequentist properties on the restricted 
parameter space of Bayesian procedures associated with tt or, more generally, truncations of the 
right Haar invariant prior measure has recently surfaced in interval estimation problems (Zhang 
and Woodroofe, 2003; Marchand and Strawderman, 2006; Marchand et al., 2008). 

In Section 2, we further describe features of the underlying model and present various expressions, 
properties, and illustrations relative to the Bayes estimators 8 n . Namely, we establish a robustness 
property, applicable to scale invariant LP loss with pit) = \t\ p , p > 0, and asymmetrized versions as 
given in ( |T3l) . stating that the Bayes estimator 5 ni does not depend on the underlying / in (J2j). 

The developments of Section 3 make use of Kubokawa's (1994) IERD (Integral Expression of Risk 
Difference) technique to derive classes of dominating (minimax) estimators of Sq(X, S) = X + CoS. 
With further analyses, which bring into play novel technical arguments of interest on their own, 
we provide several instances of (/, p) where these classes of minimax estimators include Bayesian 
estimators of the type 5 nr Namely, we establish in Sections 4 and 5 that: 

(A) The Bayes estimators 5 m with I > dominate 8q for normal models in ([1]) and general convex 
p's such that p is even. The estimator 8 no also dominates <5 for asymmetric p's such that 
\p'(u)\ > \p'(-u)\ for all u > 0; 

(B) The Bayes estimators 5 ni with I > dominate 5q for all (fixed) / in ([2]) satisfying assumption 
([3]), and whenever the loss is scale invariant L p , p > 1. The estimator 5 no also dominates Sq 
for asymmetrized versions as given in (ITBl) (where > \p(~ u )\ fo r all w > as in (A); 

(C) The Bayes estimator 5 no dominates 5 for all (fixed) / in (j2J) satisfying assumption (j3J), and 
whenever the loss is scale invariant LP with p e (0, 1). 
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The ensemble of results provide extensions of Kubokawa's normal case, scale invariant squared error 
loss result applicable to 5 m in three directions: choice of /, choice of p, and applicability to other 
Bayesian estimators 6 ni 'a. Moreover, the developments relative to (A), (B), and (C) are unified and 
contain two alternative proofs replicating Kubokawa's result. It is also notable that (C) involves 
the case of a concave in |^^| (and hence non-convex) loss. Finally, various other observations, 
including non-minimaxity results, are also given throughout the exposition and in Section 6. 

2. Preliminary results and properties of the estimator 8 n 
2.1. The underlying model 

In (j2J) and <^fy, the density of (X, U) is unimodal with central location parameter (/i, 0, . . . , 0) and 
scale parameter a. Our parameter of interest is the nonnegative fi, or median, of X, while U is 
a residual vector. Condition ([3]) is equivalent to an increasing monotone likelihood ratio (mlr) in 
(X — fi) 2 + || U || 2 of the family of densities in (j2J) when viewed as a scale family (parameter a) 
with known /x. Assumption ([3]) is, for unimodal and symmetric densities, weaker than both (a) the 
logconcavity of f(y) and (b) the logconcavity of f(y 2 ) for y > 0, with (a) implying (b), and with 
(b) equivalent to an increasing mlr property in X of the family of densities in (T5]) when viewed as 
a location family (parameter p) with known a. 

The most important and best known case covered by (j2J) and ([3]) is the normal case where (X, U) ~ 
X n+ i((/i, 0, . . . , 0), a 2 I n+ i) and f(t) oc e~*/ 2 . However, our inference results will also apply to many 
other models such as (i) exponential power densities with f(t) oc e~ atP , p > 0, a > 0, including 
Laplace densities arising for p = 1/2; (ii) the Kotz distribution with f(t) oc t m e~ at , m G (—1/2,0), 
a > 0; as well as for (iii) Student densities with f(t) oc (l + t/z/) _ ^ + ™ +1 ^ 2 , v > 1 degrees of freedom. 
The Student example illustrates a non-logconcave / (in fact, it is logconvex) which satisfies the 
weaker assumptions required here. The Student distributions, which are scale mixtures of normals, 
often serve as useful, alternative models to the normal model. Here is an interesting general situation 
for which scale mixtures inherit assumption 

Lemma 1. A scale mixture of the form f(t) = J °° vfo(tv)h(v)dv satisfies assumption (0) as soon 
as both f = fo and f — h satisfy assumption |3|). 

Proof. See Appendix. 

Remark 1. In the Student case above, both fo (a normal density) and h (a gamma density), are 
logconcave and satisfy |3|). 

Remark 2. We note that model (U|) arises for observables Y\, . . . , Y n+ i having joint density 

through an orthogonal transformation 

(Y 1 ,...,Y n+1 ) -> (X = VfiY,U 1 ,...,U n ), 

with /j = y/n0. 
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For model (j2J), (X, S = \\U\\) is a sufficient statistic with joint density fx,s on M x IR + which we 
take as equal to: 

For the normal model canonical form in (CQ), we will write the joint density of (X, S) in (EJ) as 
£0(2=<i) /»(£), with 

0(u) = (2 7 r)- 1 /2 e - 2 /2 j and ^) = j""lfoI/f-i • (8) 



r(n/2) 2"/2 



2.2. Properties of the Bayes estimators (5 



We proceed with various preliminary results, observations, and illustrations concerning the general- 
ized Bayes estimators 5 n (X, S) with respect to the improper priors in (JB}. As previously mentioned, 
one can verify that the lower bound on the power I in (jH]) guarantees that the posterior density 
of (n, a) is well defined given that (J7J) is a density. Even with a well defined posterior density, we 
further assume, and not necessarly emphasize (mainly in Sections 3,4, and 5), that the pair (/, p) 
leads to the existence of the Bayes estimator 8 n[ . 

We define form>0,«)6l,z6R, 

B m (w,z)= / / p'(u + c v + zv)v m f(u 2 + v 2 )dudv, (9) 

JO J-oo 

provided it exists. The function B m (w, z), as well as some of its properties (see for instance Lemma 
H]) will play a key role below, namely in the following representation of the Bayes estimator 8^ (X, S). 

Lemma 2. Under model §7ty, provided existence of the Bayes estimator 8 n , we have 8 n (X,S) = 
X + c S + g n (¥r)S, where g n (v) satisfies, for all y 6 R, I > —{n — \), 

Bn+i(y,g n (y)) = o. (10) 



Proof. Writing an estimator as X + c S + g(X,S), we have that the Bayes estimate 8 n (x,s) 
minimizes in g(x, s) the expected posterior loss: 



x + c s + g(x, s) — /i,, , . , 

E[p{ )\(X, S) = {x,s)\, 

a 



or, equivalently, 



00 f°° x + c s + g{x,s)-fi s™- 1 Ax-iif + s 2 1 

With the change of variables (//, a) —>-(« = ( x ~iv ; v = »^ the Bayes estimate S ni (x,s) is seen to 
minimize in g(x, s): 



J -oo 



1) 

p(u + c v + -g{x, s))f(u 2 + v 2 ) v n+l ~ x dudv . 
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Now, observe that -g(x, s) depends on (x, s) only through the function y = x/s, which implies that 
the estimator 5 n (X, S) is of the form X + c S + g ni (^r) S with g n (y) minimizing in g(y) the quantity 

p(u + c v + g(y)v)f(u 2 + v 2 )v n+l - l dudv. (11) 



Finally, the result is obtained by differentiation. □ 

We point out that g ni {y) is uniquely determined (Lemma Hj), and is a continuous function of y such 
that 

9n(y) >-V-Co for all yE R. (12) 
This must indeed be the case as the Bayes estimates S ni (x,s) are necessarily nonnegative, and 
-S ni (x,s) > •<=>■ | + Co + S , 7 r,(-) > 0. We pursue with an intriguing robustness property, 
and alternative representation, of the Bayes estimators S 7T[ for scale invariant LP loss, and their 
asymmetrized versions given by 

Pc u c 2 {t) = Ci \t\ P l(_oo,0) W + c 2 \t\ p l[o,oo)(t) , (13) 

with p > 0, ci > 0, and c 2 > 0. 

Lemma 3. For losses p Cl ,c 2 as ^ n (E2j] ; the Bayes estimators 8 n , given in LemmalE, do not depend 
on the underlying model density f provided they exist. 



Proof. From (TTTT) . we have 



oo p Vy 

2 ^ dudv 



poo p V y 

Co + 9 n (y) = argmin^ / / p CliC2 (it + hv) f{u 2 + v 2 ) v 

JO J -oo 
poo rVy u 

= argmim, / / p Cl C2 (- + h) f(u 2 + v 2 ) v^+p- 1 dudv 

Jo J-oo ' v 



r Pc 1>c2 (t+h) 

= argmm,, / — „. , — — — yttt at , 

J (1 + t 2 )( n +'+p+ 1 )/ 2 

by making use of the homogeneity of p cl]C2 and the change of variables (w,t>) — » (£ = u/v,x = u 2 +v 2 ). 
Finally, expression (TH|) tells us that S n (x, s) = x + s(co + g n (x/s)) is independent of /. □ 

This type of property seems to have first been noticed by Maruyama (see Maruyama, 2003; Maruyama 
and Strawderman, 2005) in a multivariate setting under L 2 loss. 

Remark 3. (Minimum risk equivariant estimator) 

(a) Proceeding as in the proof of Lemma\^ we obtain the useful representation X + cq(ti)S for 
the MRE estimator, with the defining equation 

OO POO 

/ p'{u + c {m)v)v m f{u 2 + v 2 )dudv = (15) 

J -oo 

for co(m), m > 1. 
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(b) A robustness property similar to Lemma\^ (also illustrated in Example^ part (C) is shared 
by the MRE estimators with respect to losses p cl , C2 an d can be established by expanding |3P 
showing that 

. [™ p Cl , C2 (t + c) 
c = argmin c J ^ ^ - t2)(w+p+1)/2 dt ■ (16) 

Example 1. (scale invariant L 2 loss, scale invariant L 1 loss and their asymmetrized versions) 

(A) For sca/e invariant squared error loss with p(t) = t 2 in Op, i/ie MRE estimator is Sq(X) = 
X , provided the second moment of X under |J|) exists. Lemma [H as we// as flTjU provide 
representations X + g n (^-)S for the Bayes estimator 5 n {X, S); I > —{n — 1). Differentiating 
|77D with respect to h, we obtain directly for y G R 

g n {y) = -E[T\T<y], (17) 

where T has density on R proportional to (1 + i 2 )-( ri + i + 3 )/ 2 . i/ere iae distribution of T is a 
multiple of a Student distribution with n + I + 1 degrees of freedom. Equivalently from [TO]) , 
we have 



oo pvy 

n+l f( n ,2 I „,2> 



B n +i(y,9n(y)) = 

M + g m {v) v ) vH+l + v ' 2 ) dudv = 

roc rvy „ f , 2 2) rf rf 

= " f™f:i v n +^f{u 2 + v 2 )dudv ' 1 ' 

illustrating the fact that the distribution of T arises as the (independent of f) distribution 
of the ratio y, with (U, V) having joint density on R x R + proportional to v n+l+1 f(u 2 + v 2 ). 
From representation (T7\ ), observe that g ni (-) decreases on R with lim^oo g- Kl (y) = (since 
uf{u 2 + v 2 )du = for all v > 0), and hence that g n (-) is positive, i.e., 5 ni expands 
on the MRE 5q. Such properties are of interest as they indicate that the amplitude of the 
expansion 5 ni (x, s) — 8o(x, s) decreases in x for fixed s, and increases in s for fixed x (in fact 
(S ni (x,s) — 5o(x,s))/s increases in s). Such a property resonates back to Katz (1961) where 
in the normal case with known a, the Bayes estimator with respect to a flat prior for p, on 
(0, oo) expands X by the amount a which decreases in x and increases in a. Below, 
we establish such properties for general convex p in Lemma 0, as well as scale invariant LP 
concave loss with p G (0,1) in Lemma{5[ Finally, we point out that alternative expressions 
for 5 ni in the above normal case were given by Kubokawa (2004), as we ^ as Marchand, Jafari 
Jozani, and Tripathi (2011). 

(B) As above, for scale invariant absolute value error loss with p(t) = \t\ in Op ; the MRE estimator 
is S (X) = X . For I > —{n — \), S n (X, S) = X + g n (^)S is obtainable from |7jj ) yielding 

-1 fFn+i(y). 
2 

where F m and F,~ x are the cdf and inverse cdf of T having density on R proportional to 
(1 + t 2 )~( m + 2 )/ 2 . As above, it is easily seen directly that such a g m {') decreases on R, that 
lim^oo g m (y) = 0, that S m (x,s) expands once again on 6 (x,s) for all (x,s) G R x R + , and 
the difference between these estimates decreases in x/s. 



g ni (y) = -median\T\T < y] = -^S^l) , (19) 
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(C) Consider now asymmetrized L l losses p Cl ,c 2 i n |Z2P with p = 1. 5?/ making use of Remark 
0, £/ie MRE estimator is given by S (X) = X + c S, with c independent of f , and c (n) = 
— F~ 1 ( c ^ ) and F" 1 £/ie inverse cdf given in part (B). Fori > — (n — 1), we obtain from (Lfy 
S n (X, S) = X + c (n)S + g n (§)S with g n (y) = -c (n) - F^(^^F n+z (y)) ; i/ras extending 
/ flgj) which occurs for c\ = c%. Observe here that lim^oo g ni (y) = —co(n) + co(n + I), which 
does not equal in general, the exception being precisely 1 = 0, and/or c\ = C2- This property 
is more general as seen below in LemmalB 

We pursue with further properties relative to B m (-,-) and g Wl (applicable when these quantities 
exist). 

Lemma 4. For all a > 0, y G M, / > —{n — 1), and strictly bowled-shaped p, 

(a) B n+l (y + a,g n (y)) > 0; 

(b) B n+ i(y,z) is nondecreasing in z whenever p is also convex; 

(c) Hindoo B n+ i(y, 0) = whenever I = 0; or whenever I ^ and p is an even function. 

Proof. Part (b) is obvious given the convexity of p, while part (c) follows from the given repre- 
sentations ( !T5|) and ( TTUj) . For establishing (a), suppose, in order to arrive at a contradiction, that 
B n+ i(y + a,g ni (y)) < 0. This would imply C\ < 0, where 

Ci = / / p'{u + c v + g Vl (y)v) v n+l f(u 2 + v 2 ) du dv . 

JO Jvy 

Now, observe that for (u,v) G I{u,v) = {(u, v) : vy < u < v(y + a)}, we have by ( IT2l) : u + cqv + 
g n (y)v > vy + c v + g n (y)v > 0, implying p'(u + c v + g n (y)v) > 0, (for such (u,v)'s G I(u,v)). 
This renders C\ < impossible, and yields the result. □ 

The strictly decreasing property of g m that follows in Lemma [5] is a critical property that we will 
exploit later for the risk comparisons. We do not know how far the property can be extended 
for non-convex p, but we do establish here, and use later, such a property for L p losses and their 
asymmetrized versions for the non-convex choices p G (0, 1). 

Lemma 5. For I > — (n — 1), 

(a) g %l (y) is strictly decreasing in y whenever p is convex; 

(b) g n {y) is strictly decreasing in y whenever the loss is p Cl ,c 2 as ^ n (123) with p G (0, 1). 

(c) For strictly bowled-shaped p, lim^oo g n {y) = —co(n) + co(n + I), where co(m) is defined in 
173]) . Consequently, lim^oo g n (y) = whenever 1 = 0, or I ^ and p is even. 
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Proof, (a) It suffices to show that we cannot have g^Xv + e ) — d^iiv) f° r some y G R, e > 0. 
Indeed, if this were the case, it would follow, using defining equation ffTUl) and part (a) of Lemma 
H that 

= B n+i (y + e, g ni (y + e)) > B n+l (y + e, g^y)) > , 

which is not possible. 

(b) Set s(y) = — Co — g^iy) and rewrite representation (p3J as 

s(y) = argmin s £[p Cl , C2 (T - s)\T < y], (20) 

with T having density proportional to (1 J r t 2 ^- ( J l + l +p+ 1 )/ 2 Q n R. Observe that the family of densities 
for T\T < y has strictly increasing monotone likelihood ratio in T with parameter y. Now, consider, 
for ai < a 2 , the function p Cl , C2 (t ~~ a i) ~ Pci,c 2 (^ — a 2)? which changes signs once from — to + as a 
function of i as t increases on R, and infer that 

H(ai,a 2 ,y) = E[p CUC2 (T - ax) - p cl , C2 (T - a 2 )} 

has a single root, and changes signs once from — to +, as a function of y, as y increases on M, 
given the mlr property (e.g., Lehmann, 1986). Suppose now, in order to arrive at a contradiction 
that g n is not strictly decreasing, i.e., s is not strictly increasing and there exists y% < y\ such 
that a.2 = s{y 2 ) > s(yi) = a\. Then, we would have with the definition of s(y) in f l20|) and the 
properties of H: H(s(y\), s(y 2)^2) > and H(s(yi),s(y2),yi) < which leads to a contradiction 
and establishes the result. 

(c) This follows by matching expression ffTU]) when y — y 00 with (|15|) . □ 

Remark 4. The above proof in (b) goes through for all losses p Cl)C2 > including the convex cases with 
V > 1- 

The following results permit the ordering of Bayes estimators S ni in terms of the power I in the prior 
measure 717 in (jSJ). 

Lemma 6. For the normal model in y G R, and convex and even p, the quantities g n {y) 
decreases in I, I > — (n — 1), provided they exist. 

Proof. See Appendix. 

Corollary 1. For models (0) / satisfying assumption §3±), t/Gl, and sca/e invariant LP loss 
with p > 0, g n (y) decreases in I, I > —{n — 1), provided existence. 

Proof. Lemma |3] tells us that g n (y) is independent of / and thus matches the normal model g n {y) 
and Lemma [6] tells us that such g ni (yys decrease in I whenever p is even as for the L p loss here. □ 

3. Minimax Conditions for general p and / 

For estimating p > in (JSJ) or in (JTJ) with unknown a > under strictly bowled-shaped loss p(^—^), 
we establish here useful sufficient conditions for an estimator 5(X, S) to be minimax. We first 
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make use of Kubokawa's IERD technique in Theorem [TJ Proposition [T] (below) then extracts a sign 
varying condition for minimaxity which will serve as the basis for further analysis for the specific 
cases of normal models and general convex p in Section 4, and for LP losses and their asymmetric 
versions p Cl , C2 with general / satisfying assumptions (J3J) in Section 5. Various other technical results 
and remarks, including a condition for non-minimaxity with applications, are also introduced in 
this section. We consider the following subclass of scale invariant estimators. 

Definition 1. C = {5 g (X,S) : S g (X,S) = 5o(X,S) +g(^)S, with g absolutely continuous a.e., 
nonincreasing , non- constant, and lim^oo g(t) = 0}. 



These estimators in C expand upon 8q, in view of the restriction p > 0, include 5 W0 and the 
generalized Bayes estimators 8 m ] I ^ 0, I > — {n — 1); for even p as seen by the properties given 
in Lemma Under invariant losses as in (J3J), such estimators will have frequentist risk R(6,S g ) 
depending on 9 = {p,cr) only through the maximal invariant A = p/o~, and we seek conditions for 
which such a risk falls below the constant risk of the MRE estimator So for all A > 0. As mentioned 
above, such improvements will necessarily be minimax estimators since So is minimax. Hereafter, we 
will just refer, for the most part, to such improvements as being minimax estimators. The focus is 
largely on the generalized Bayes estimator S no , which will be seen to be minimax for various settings 
of (/, p) and which provides a benchmark in the sense that estimators S g G C will be minimax for 
convex p under the simple condition that 5 g not expand on So as much as (Theorem [TJ (ii)). 
In turn, for various choices of (/, p) with p even, and by appealing to Lemma El these classes of 
minimax estimators will contain the generalized Bayes estimators 5 ni 's, I > 0. We now pursue with 
an intermediate dominance condition. 

Theorem 1. For estimating p in (TJ|) or (Q) with p > 0,cr > ; an estimator 5 g G C is minimax, 
under strictly bowled shaped loss p(^-^) whenever either one of the following conditions holds for 
all A > and y G {y : g'(y) < 0}: 

(i) 

oo pvy—\ 

/ p'{u + cqV + g(y) v) v n f(u 2 + v 2 ) du dv < , 

J-oo 

or 

(ii) p is convex, g < g no and ip p (\,y) < 0, where 

poc pvy—X 

ip P (X,y)= / p'{u + Cqv + g no (y)v) v n f(u 2 + v 2 ) du dv . 

JO J -oo 

Proof. With p'(-) increasing by the assumption of convexity, condition (ii) implies (i) so that we 
only need to establish the sufficiency of (i). Following Kubokawa (1994), write for S g (X, S) G C, 

x + c s-p x + c s + g(£)s - p x + cos + g(y)s - p 
P( ) - P{ ) = P( )\ y y=x/s 

s , x + c s + g(y)s- p , 

-p ( )g {y)dy ■ 

x/s ° ° 
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Now, use the above expression for the difference in losses to write the difference in risks at 6 = (//, a) 
as: 



= R(6,5 )-R(6,5 g ) 



a 



/ s i 9{y)p{ )dy}f x ,s{x,s)dxds 

J0 J-oo Jx/s ° 

9{y){ / P{ )-^f( 3 )dxds}dy, 



{g'(y)<0} JO 



(21) 



since g' < a.e. Now, the difference in risks A g (9) will be nonnegative for all 6 G as long as for 
all y G E such that g'(y) < 0, fi > 0, a > 0, the bracketed term in (121]) is less than or equal than 0, 
which is equivalent to (i) with the change of variables (x, s) — > (u — ^^,v = -). □ 

Remark 5. Notice that ip p (0,y) = B n (y, g„ (y)) = for all y G R by virtue of the definition of 
9w in (G2P- Therefore, the risks of S no and 5 match at the boundary of where /i = 0, a > 0. 
Moreover, if S g expands more that 8 no ( whether or not 5 g G C), then the risk at the boundary of 5 g 
will exceed that of 5q, hence giving a condition for non-minimaxity. This is so given that 

R((0,a),6 g ) = E m (p(6 g (X,S))) 
>E m (p(5 no (X,S))) 
= R((0,a),5 7TO ) 
= R((0,a),S ), 

since S no (X,S) > with probability one, and p is increasing on (0,oo). As a consequence of the 
above, and of Lemma and Corollary Ui we have the following non-minimaxity result. 

Corollary 2. For estimating p in (TJ|) or (RD with p > 0, o > 0, the generalized Bayes estimators 
5 ni with — (n — 1) < / < are not minimax whenever (a) / is normal and p is even and convex, or 
whenever (b) / satisfies assumption (dj] and the loss is invariant LP with p > 0. 



Analogously, we point out that 8^ does not dominate any other minimax estimator 8 g G C taking 
nonnegative values and satisfying (ii) of Theorem^ since such S g 's shrink 5 no and R((Q,a),8 g ) = 
E m) (p(8 g (X,S))) < E m (p(8„ (X,S))) =R((0,a),8 nQ ). 

Remark 6. A plausible alternative to the MRE estimator 8$ is, of course, its truncation %(X, S) = 
max(0, 8q(X, S)). Clearly 8fi improves upon 8o for bowl shaped p, since for all p > 0, a > 0, 

(x, s) G 1 x R + , with strict inequality occurring with positive 
probability. Moreover, the estimator 8$ belongs to the class C with g${y) = max(0, — y — cq), and 
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satisfies condition (i) of Theorem\J\ with {y : (go)'(y) < 0} = (—00, — cq) since 

foo pvy—X 

p'(u + cqv + g^{y) v) v n f{u 2 + v 2 ) du dv 



j —00 

00 pvy—\ 

\ p'{u - vy) v n f{u 2 + v 2 ) du dv 
'0 J -00 

/*oo rvy—\ 

< / p'(-A) r/V(w 2 + r 2 ) (iuA 
<0, 

/or a// A > 0. Finally, the observations of Remark^ apply to 8q, with 5q a shrinker of 8^, and S no 
not dominating 5$. 

With Theorem^ our attention focuses on the quantity ip p (\, y) and testing the condition i/j p (-, •) < 
on R + x R for various choices of p. Now, since 

^ (0, y) = 0, and lim ip p (X, y) = for all y G R, (22) 

A— >oo 

VV>('>2/) cannot be monotone on [0, 00) for any p and y G R. We are thus led to analyzing the 
behaviour of J^ijj p (\,y). 

Proposition 1. Let k(y) = y + c + gir {y), fx,y(t) be a Lebesgue density on (0, 00) proportional to 



t U f ( A 2 (l + y 2 ) {(t - -^) 2 + V7~~" 2\2 ^ 

V i + r (1 + r) 2 

l/fc(y) 



and 

D p (A,y) = / \ P '(X(tk(y)-l))\f x , y (t)dt- I \p'(X(tk(y)-l))\f x , y (t)dt. (23) 

JO Jl/k{y) 

Suppose further that D p (-,y) changes signs once from — to + on [0, 00) for all y G R. Taen, /or 
estimating p in (TJ|) or (?p under assumption (|3j] p > 0, a > ; 

(i) £ae generalized Bayes estimator 5^ is minimax, for strictly bowled shaped loss as Zona 
as o'tto G C; 

(ii) /or o~ 9 G C, £ae condition g < g^ is sufficient for 5 g to be minimax under convex loss p(^p). 
Proof. We have 

— i> p {\ V) = -J P'(vk(y) - A) v n f((vy - A) 2 + v 2 ) dv 

POO 

oc -A" +1 / p'(\tk(y) - A) f X;y (t) dt (24) 
Jo 

<xD p (\,y). (25) 

Therefore, under the given assumptions on the sign changes of D p (-, y), we infer that, for all y G R, 
ijj p (\,y) decreases, then increases as A varies on [0, 00). Finally, the result follows from Theorem [T] 
and property ( |22l) . □ 
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Remark 7. (i) From (2$) , note that 

q^ p {\ 2/)U=o+ = -] p'(vk(y)) v n f(v 2 (y 2 + 1)) dv < 0, 

since k(-) > from ( TIP)) , and p'(-) > on [0, oo). Hence, Proposition U\'s sign change 
assumption on D p (-,y) is consistent, for any strictly bowled-shaped p, with the behaviour of 
ip p (\, y) for X near 0. 

(ii) Turning to the families of densities {f\, y {-), A G [0, oo), y G K} ; they can be shown for y < to 
possess a decreasing monotone likelihood ratio (mlr) in T, or equivalently inW = (T — t^i) 2 > 
with A viewed as the parameter. Indeed, for Ai > A > ; setting c*i = X 2 (y 2 + 1) and 
e — (y 2 + I) -2 ; we have 

h ,y(t) f{a {w + e)) 

2 

which decreases in w, w > t^^jw, given assumption 

(iii) In the normal case, the densities fx, y (-) may described as weighted (by the factor t n ) positively 
truncated N{y/{1 + y 2 ), 1/(A 2 (1 + y 2 ))) densities. They have been recently studied in related 
work of Marchand, Jafari Jozani, and Tripathi (2011) where quantiles are estimated under 
the restriction p > 0. 

We conclude this section with a very useful technical result. 

Lemma 7. Let k(y) = y + Cq + g-K {y) as in Proposition^ and let p be either an even function, or 
more generally satisfy \p'(— u)\ < \p'(u)\ for all u > 0. Then we have > max{0, x+p-}- 

Proof. The positivity of k(y),y G M. follows from ( TT2|) . To establish that > j^s, we assume 
the contrary and show that this would imply D p (X,y) < for all A > which is not possible given 
( 12^1) and ( 1251) . Indeed, we would have, for all A > 0,y G M, under the given assumption on p 

D P (X,y) 

r l/k{y) p2/k(y) i 

< / \ p '(Xk(y)(t-—))\f x , y (t)dt- \ p >(Xk(y)(t-—))\f Xjy (t)dt 
Jo k(y) Ji/k{y) k{y) 

\p'(Xk(y)(t - _))| (f x , y (t) - f x>y ( _ - t)) dt 

<o, 

given that f X , y {t) < f X , y {^ - 1) for all t G (0, l/k(y)) whenever < j^. □ 

The inequality > max{0, j^p} will be exploited as a technical result, but it also provides an 
interesting upper bound for the generalized Bayes estimator 8 no , namely 

2 

X S 

S^tx, s) = s k(-) < x H , for x > 0, 

s x 

applicable to all pairs (/, p) for which 5 no exists, with / satisfying ()3]), p satisfying the conditions 
of Lemma [71 
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4. Minimax results for the normal case 



Here is a minimax result applicable in the normal case, to Bayes estimators 5^, and for general 
convex losses that are either even functions or, more generally, that penalize the rate of over- 
estimation more sharply than the rate of underestimation in the sense 

\p'{-u)\ < p\u) , for all u > . (26) 

Theorem 2. For estimating p, in the normal case in with p > 0, a > under convex p in 



(a) the condition g < g no suffices for an estimator 5 g e C to be minimax in cases where p satisfies 
condition 



(b) such minimax estimators include the generalized Bayes estimator 8^ under losses p satisfying 
l[2b}) . and all S Wl with I > when p is even. 



Proof. Given that 5 ni € C for I = 0, and for I > when p is even by virtue of Lemma [5], the first 
part of (b) is simply a restatement of (a) for the Bayes estimator 5 W0 , while the part relating to 8 n 
with I > follows also from (a) and Lemma O The rest of the proof concerns part (a) and we apply 
Proposition [H From (|23l) . we have with the change of variables u = X(tk(y) — 1): 

D ^-jm E[ - p ' m ' 

where U has density proportional to 

^<5 + 1 >*y i <-W«>- (») 

Since —p' changes signs once from + to — on R, a decreasing in u monotone likelihood ratio property 
of the densities in (|27|) with respect to the parameter A will suffice to establish that D p (\, y) changes 
signs from — to + on [0, oo) as a function of A > and permit us to apply Proposition Now, 
the densities in (j27|) may be written as 

h X { U ) OC (^)" / (C(U + A&) 2 + d\ 2 ) 1(_a,oo)(w) , 

with c = (1 + y 2 )/k 2 (y), 6=1- (yk(y)/(l + y 2 )), and d = (1 + y 2 )~ l . Notice that we have c > 
by virtue of (fl2|) . and b > by assumption fl26|) and Lemma [71 Finally, in the normal case with 

f(t) = (27r) (n+1)/2 e"' /2 , tha ratio is, for A x > A > 0, undetermined for u < -Ai, equal to 

+oo for k G (— Ai, — Aq], and otherwise proportional to 



u + Ai 

M + An 



— 6«i(Ai— Ao) 



2 It is interesting to point out that the arguments here apply as well to strictly-bowled shaped losses. As well, 
only a stochastic increasing property for the densities / is required. The monotonicity of however is guaranteed 
by the convexity of p (Lemma [SJ , which is assumed here. 
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which is indeed decreasing in u for u > — Aq, and which establishes the result. 



□ 



The normal case minimax results of Theorem [2] in part (a), and applicable to the generalized 
Bayes estimator 5 no , were previously obtained for the specific case of scale invariant L 2 loss by 
Kubokawa (2004). He works directly with the Bayes estimator in Example 1 to derive the key 
required analytical properties, namely the monotonicity of g %0 in Lemma [5] and inequality (i) of 
Theorem [TJ With Kubokawa's analysis specific to scale invariant L 2 loss, the normal model and 
the estimator 5 no , our unified development above contrasts and provides extensions with respect to 
the loss and the prior. In the next section, we give extensions with respect to the model for scale 
invariant LP losses and asymmetric versions. 

5. Minimax results for scale invariant L p losses and their 
asymmetric versions 

The minimax results of this section are applicable for the wider class of models, or choices of /, in ([2]) 
with assumptions ([3]). As well, these findings concern scale invariant LP losses p > 0, and the 

more general p cltC2 in CE2D with c 2 > c\. For these losses, (j2"3"j) reduces to D p (X, y) = p\ p ~ x Ex[g y (T) ], 



with T ~ f X) y and g y (t) = c x (1 - tk{y)) p 1 l(o,i/fc(»)(t) - c 2 (tk(y) - l) p 1 l(v*( v ), 00 )(*)- With this 



representation, observe that g y (-) changes sign once on (0, oo) from + to — , so that E\[g y (T)] 
changes signs from — to + as A varies on [0, oo), in view of sign change properties and the mlr 
property of Remark [7] (ii). Therefore, D p (X,y) varies indeed, as a function of A G [0, oo) from — 
to + as prescribed in Proposition [TJ for y < and losses p cl)C2 - For y > however, the situation is 
more delicate. We continue with the non-convex case with p e (0, 1), and this will be followed by 
the convex case with p > 1 . 

Theorem 3. For estimating \i in (TJ|) or |7j) with p, > 0, a > under scale invariant LP loss \^^\ p 
with p G (0, 1), the generalized Bayes estimator 5 no is minimax. 

Proof. With 5 no G C by virtue of Lemma we seek to apply part (i) of Proposition [1] to show 
that 5^ is minimax. For p(t) = \t\ p with p > 0, we reexpress f l23|) as 



D P (X, y) oc E[A(T) B(T)] = E[G(S)] 



with 



G(s) = E[A(T) B(T)\S = s), A(t) = t n \t 



1 



p-i 



B{t) = -1 + 21 (0 i ,(*), 

V <k(y)> 



k(y) 




and 



S = (T 



y 




l 



l + y 



(i + y 2 ) 2 ' 
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Given assumption ([3]), the family of densities of S are seen to have a decreasing monotone likelihood 
ratio in S with parameter X 2 (l + y 2 ). So, in accordance with Karlin's sign change analysis, to prove 
the result, it will suffice to show that 

G(s) changes signs once from + to — as s varies on (-. — , oo) (28) 

(i + zrr 

to establish that D p (X,y) changes signs as prescribed by Proposition [TJ We proceed by treating 
separately the cases: (i) — > j^z an d (ii) < — < j^p- Here, we have made 
use of Lemma [7J to discount the remaining possibility — < 0. 

Case (i): Set s = - j^) 2 + Observe that, whenever s > s Q , P{T > -^\S = s) = 1 

implying P(B(T) = -1\S = s) = 1 and G(s) < 0. Similarly, if s < s , then P(B(T) = 1\S = 
s) — 1 and G(s) > 0. Hence, the above establishes (128]) for case (i). 

Case (ii): Here, we set Si = nj^p + j+^si so that s < Si. As in (i), we verify that G(s) > for s < s , 
and G(s) < for s > s\. Finally, for s G (so,si), the conditional distribution of T\S = s is 
a two-point uniform discrete distribution on {ti, ^2}, with t\ = + A, t 2 = — A and 



A 



1 



(l+J/2)2 

We hence obtain 



G{s)= l -{A(t 2 )-A(t 1 )) 



+n I. ip-1 ,n I . 



< 



since t 2 < h, n > 1; |*i - ^| = - ^ + A < + ^ + A = \t 2 - and p - 1< 0. 

Hence, the above establishes ( )28l) for case (ii) and completes the proof. □ 

Theorem 4. For estimating p in (Ej) or |?|) under assumptions (TJ}) ; u>i£/i /i > 0, a > and with 
loss p Cl ,c 2 > P > 1 an d c 2 > ci, 

(a) £/ie condition g < g no suffices for an estimator 5 g G C to be minimax; 

(b) such minimax estimators include the generalized Bayes estimator 8 no , as well as all generalized 
Bayes estimators 5 n , I > for the symmetric case c\ = c 2 . 

Proof. For losses p as in (fT3|) . we may write 

fkb) f°° 

D p (X,y)=p\ p {d h x , y {w)dw-c 2 h x , y {w)dw}, 

k(y) 

with h\ y {-) a probability density function on (0, 00) proportional to \wk(y) — l| p_1 f\ y {w). From 
this, we see that D p (X,y) is positive iff P\(W > l/k(y)) < c\/(c\ + c 2 ), where W is a random 
variable with pdf h\ y . We show below in Section [7^41 of the Appendix that, whenever > j^z, 
the quantity P\(W > l/k(y)) decreases in A on [0,oo), which means that D p (-,y) changes signs 
from — to + on [0, 00). The result then follows from Proposition [T] and Lemma UJ □ 
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6. Concluding Remarks 



We have considered the problem of estimating a lower bounded location parameter for a wide array 
of spherically symmetric location-scale models with a residual vector as represented in model (T5J), 
with unknown scale, and under scale invariant loss as given by (jl]). With a relative paucity of findings 
for such problems when the scale parameter is unknown, we have established the minimaxity of the 
generalized Bayes estimator 5 no for normal models and convex loss, as well for more general models 
and scale invariant L p loss and asymmetric versions p Cl , C2 given in ([TBI . Moreover, we have shown 
the role of 5 no to be pivotal, in the sense that it provides an upper threshold condition necessary 
for the minimaxity of many estimators. Other minimax estimators are also obtained, including 
generalized Bayes estimators 5 1Tl when I > and the loss is convex and even in the above situations. 
The results represent extensions of Kubokawa's results (2004) applicable to scale invariant L 2 loss. 
Much of the treatment is unified and exploits general features of the model and the loss with incisive 
analysis and novel representations. Various other observations are given, including the robustness 
of the Bayes estimator 5 7Tl with respect to the choice of / in model fl2]). 

As illustrated by Marchand, Jafari Jozani and Tripathi (2011), the normal case improvements pro- 
vided for scale invariant L 2 loss yield applications for two-sample problems where Y$ ~ N(fii, a 2 ); i = 
1,2 with unknown /i l5 /i 2 , a 2 , where the objective is to estimate /ii (or fi 2 ) with the additional in- 
formation of the ordering [i\ < /i 2 . Despite the advances presented here, minimax extensions to 
other strictly bowled-shaped losses, although plausible, are still lacking. Furthermore, numerous 
questions remain unanswered such as the admissibility of the above minimax estimators, the in- 
vestigation of wider classes of Bayes estimators for minimaxity, and related tests of minimaxity for 
multivariate location-scale problems with order restrictions. 

7. Appendix 

7.1. Proof of Lemma [I] 

We have 



where Z has density proportional to fo(z)h(z/ 1) on IR + . Now, observe that this scale family of 
densities for Z has increasing monotone likelihood ratio in Z, with parameter t, as a consequence 
of assumption ([3]) for h. Finally, the result follows from representation f|29|) with this monotone 



tf'(t) 



fit) 




J °° tv 2 f (tv) h(v) dv 
f™vf (tv) h(v) dv 

f™ z fo( z ) h (z/ 1) dz 
JcT ' fo(z)h(z/ 1) dz 



(29) 



likelihood ratio and since 



decreases in z by assumption ([3]) for f . 



□ 



M*) 
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7.2. Proof of Lemma [6] 



We fix y 6 R throughout and set Co = given that p is assumed even. First, observe that by 
differentiating (fTUl) for the normal case with f(u 2 +v 2 ) = 4>{u)h(v) in (jSJ), we have -^jB n+ i(y, g n (y) = 
which implies 

d { r vy 

p' '(it + g n (y)v)<f)(u)du 



dl 

vy 1 

p'(u + g n {y)v)(j)(u)du > \og(v)]v n+l h(v)dv = 0. 

-oo J 

Hence, given that p 1 is increasing, to show that g ni {y) decreases in I, it will suffice to show that 
/ > 0, where 

poo pvy 

1= (\ogv)A y (v)v n+l h(v)dv, aud ^(u) = / p'(u + g n (y)v)(/)(u)du . (30) 

JO J-oo 

Now, we will show below that 

A y (v) changes signs once as a function of v from— to + . (31) 

Applying Lemma [H which is stated in the Appendix, with £ ~ £ n+l h(£) l(o,oo)(£)> r (0 = l°g(£)5 
and s(£) = A y (£), we infer that / > since E[A y (£)] = given the definition of g^y) in ffTUj) . 
There remains to establish ( |3T1) . which we proceed to do separating the cases: (i) y < and (ii) 
y>0. 

(i) Case y < 0. Let Vq be such that A y (vo) = 0. Such a value exists since the average value of 
A y (£) under the above density for £ is equal to 0. For e > 0, we have 

A y (v + e) 



$((v + e)y) 

r{vo+e)y ^r u \ 
= J o P'{u -ey + e(y + g ni (y)) + g n (y)v ) ^- + ^ 



Jo n(^o + e)y) 

= ^(vo, e) (say) , 

with equality if and only if e = 0, given f lT2|) and since p' is increasing. Now, observe that the 
ratio of densities ^^^^/^(vo+y)^ 1S mcreasm g m u ' ' ■> f° r u ' e ( — °°; v oy) and ey < 0. Hence, 
this monotone likelihood ratio property implies that ^^f^k > C y (vo, e) > C y (vo, 0) = 0, for 
e > and y < 0, yielding (13T1) for y < 0. 



(ii) Case y > 0. As in (i), let vq be such that A y (vo) = 0. Using this, as well as property (|T2l) . the 
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nonnegativity of g n (-) (Lemma [5]) (since Co = 0), and the convexity of p, we have for e > 0: 

/(v +e)y 
p'{u + {g n (y)){v + e)}(p(u)du 
-oo 

> A y (v ) + / p'{u + (g n {y))v } (j){u)du 

J v y 

> 0. □ 



7.3. Lemma used within the proof of Lemma [6]. 

The following result is well known and its proof is left to the reader. 

Lemma 8. Let £ be a continuous random variable, letr(-) be a continuous and increasing function 
on the support of and let s(-) be a continuous function which changes signs once from — to + at so 
on the support of £. Then, we have E[r(£) s(£)] > r(so) £7[s(£)] ; and, in particular if E[s(£)] = 0, 
then E[r(£) s(f) ] > 0. 

7.4. Theorem [H proof of a monotonocity property for P\(W > l/k(y)) 
We wish to show that 

P\(W > ) decreases in a whenever -r-^-r > a, (32) 
k(y) k(y) 

with W having pdf h\ >y (w) on (0, 00) proportional to \wk(y) — l| p_1 w n f(a{(w — a) 2 + e}, a = j^t, 
a = A 2 (l + y 2 ), and e = (y 2 + 1)~ 2 . We have 

d -p,iw> 1 * 



o\ ^ k( y y 

d \wk(y) - l\ p - 1 w n f(a{(w - a) 2 + e} dw 

= 2A(1 + y 2 )— < 

act J Q |wA;(y) — l| p i w n /(a{(w-a)' i + e}aw 

<=► J3[ 7 (a{(W - a) 2 + e})|W > -j^A > EfriaKW - a) 2 + e})}, (33) 

under pdf h\ >y , with 7(2) = . Taken together, the following hence form a sufficient condition 
for (E3D to hold: 

(i)E[j(a{(W - a) 2 + e})\W > -^-} > E[j(a{{W - a) 2 + e})\W < a] , and 

k{y) 

(ii) E[-f(a{(W - a) 2 + e})\W > — T ] > E[-f(a{(W - a) 2 + e})\a < W < ' ' 



%) %) 

Condition (ii) is immediate since ^{a{{W — a) 2 + e}) increases in (W — a) 2 on (a, 00) by assumption 
(E]). For (i), set Z = \W-a\ so that Z\W > l/k(y) has pdf proportional to \(a + z)k(y) - l| p_1 (a + 
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z) n f(a(z 2 + e)) ln/k(y)-a,oo)(z), while Z\W < a has pdf proportional to |1 — (a — z)k(y)\ p 1 (a — 
z) n f(a(z 2 + e)) 1/Q >a \(z). We thus have the ratio 

oo if z < a, z < — a; 

fz\w<a{z) Jo ifz>a,z>^-a; 



/z|W> T f T ( 2; ) /a-z\ n /^2A:(?/)+(l-afc(s/))\ P 1 -r n ^ ] 

( (irMlj 1 fO<z<a,z> W) -a. 

Since both 2=5 and decrease in z for z < a and 2 > l/k(y) — a, with 6=1 — aA;(y) > 

(Lemma [7]), we have a decreasing monotone likelihood ratio. Finally, with 7(a(z 2 + e)) increasing 
in z > by (jHj), condition (i) follows and our proof of ( |32l) is complete. □ 
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